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Abstract: We predict a sharp crossover from nonlinear self-defocusing 
to discrete self-trapping of a narrow Gaussian beam with the increase of 
the refractive index contrast in a periodic photonic lattice. We demonstrate 
experimentally nonlinear discrete localization of light with defocusing 
nonlinearity by single site excitation in LiNbC>3 waveguide arrays. 
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1. Introduction 

The growing interest in the study of nonlinear periodic structures is motivated by the fact that 
novel physical phenomena can be observed due to the interplay between nonlinearity and pe- 
riodicity [1]. The presence of periodically varying refractive index in the medium results in a 
bandgap structure of the transmission spectrum. This subsequently affects the propagation of 
optical beams, which is determined by the dispersion curves of different transmission bands and 
forbidden gaps. By engineering a periodic structure, it becomes possible to manage the strength 
and type of wave diffraction [2], and therefore control the self-action of light in nonlinear me- 
dia [3]. If the refractive index decreases with light intensity due to nonlinear response of the 
material, the beam normally experiences broadening due to the self-defocusing. However, in 
periodic photonic structures, the same type of nonlinearity allows for beam localization, and 
this effect is analogous to nonlinear self-trapping in discrete systems [1]. 

In this paper, we study the self-action dynamics of an initially narrow Gaussian beam prop- 
agating in a nonlinear defocusing medium with periodically modulated refractive index. We 
demonstrate that there is a sharp transition from self-defocusing to discrete self-trapping when 
the depth of index modulation (An) is increased. We also observe experimentally the formation 
of a self-trapped state from a single- site excitation in the defocusing regime. 

2. Discrete self-trapping in photonic lattices 

Nonlinear propagation of light in photonic structures with a periodic modulation of the optical 
refractive index in one transverse spatial dimension, such as waveguide arrays [1, 3, 4, 5, 6, 7] 
or optically-induced photonic lattices [8, 9, 10, 11], is commonly described by the so-called 
tight-binding approximation [12, 13], 

j^l + p an + C(a n -i+a n +i) + Y\a n \ 2 a n = 0, (1) 
az 

where z is the propagation coordinate, a n (z) is the mode amplitude in the n— th waveguide 
[Fig. 1(a)], and /3 is the propagation constant. Coefficient C stands for the nearest-neighbor 
coupling between the waveguides, and the last term in Eq. (1) accounts for mode detuning 
through the intensity-dependent change of the refractive index. 

In the case of self-focusing nonlinearity, Eq. (1) supports the so-called unstaggered discrete 
solitons [1, 13]. These are strongly localized states where the light is confined to a few waveg- 
uides, and the amplitudes of the modes of neighboring waveguides are in-phase. Importantly, 
defocusing nonlinearity can also support discrete self-trapping of light, but in the form of stag- 
gered solitons when the amplitude of the modes in neighboring waveguides is out of phase [14]. 
This can be understood by noting that the form of Eq. (1) remains unchanged after the trans- 
formation: a n — > (—l) n a* and y — > — y. It follows that the beam dynamics in the framework of 
Eq. (1) is fully equivalent for positive (y > 0) and negative (y < 0) nonlinearities, with the only 
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Fig. 1. (a) Schematic of a discrete lattice with nearest-neighbor coupling, (b) Refractive 
index profile of coupled LiNb03 waveguides. 



difference being in the phase structure, provided only a single site is excited at the input, i.e. 
a n (z = 0) = for n ^ 0. In particular, when the input intensity is high enough, a discrete soliton 
should form for either type of nonlinear response. 

3. Nonlinear localization in continuous periodic photonic structures 

The tight-binding model (1) predicting a universal self-trapping scenario is valid when the op- 
tical field can be represented as a superposition of weakly overlapping modes of the individual 
lattice sites. Such condition is generally satisfied when the refractive index contrast in a pe- 
riodic structure is sufficiently large [15]. However, in the case of defocusing nonlinearity, the 
beam dynamics changes dramatically for smaller index contrast. Indeed, if the lattice modula- 
tion is absent or very weak, the beam will only experience nonlinearly enhanced defocusing 
instead of self-trapping. In order to study the crossover between beam self-defocusing and dis- 
crete self- trapping, we model the beam propagation using a normalized continuous nonlinear 
Schrodinger equation [16] for the slowly varying field envelope E(x,z), 

dE d 2 E 

i— + + ^(\E\ 2 )E + p An(x) E = 0, (2) 

oz ox z 

where D = z s A/(4;rnoJtJ) is the diffraction coefficient, p = 27lz s /k, x is normalized to x s , z 
is normalized to z s > We choose the parameters to match the conditions of our experiments 
described below. The linear refractive index change An is taken as An(x) = £E n exp[— (x — 
n d) 2 /w 2 ], where t; defines the modulation depth. We take the values of waveguide width 
w = 12/im and the waveguide spacing d = 19/im, and then the corresponding refractive in- 
dex contrast is An max = 0.442 We note that this index profile [Fig. 1(b)] is defined by the 
experimental realization, however we have verified that our conclusions are valid for different 
lattice profiles. Other parameters are: A = 0.532/im, no = 2.234, x s = 1/im, z s = 1mm, and 
&(t) = 1.5(1 -f for photovoltaic defocusing nonlinearity [17]. 

The periodic modulation of the refractive index results in the formation of a bandgap 
structure for the wave- vector components, as the light propagates along the waveguides. In 
Figs. 2(a,b) we present the calculated bandgap spectrum (see Ref. [18] for details of such cal- 
culations) for the waveguide array shown in Fig. 1(b) for two values of the refractive index 
contrast An max . In the case of defocusing nonlinearity, solitons can form at the bottom edge of 
the first band, and their propagation constant is shifted deeper into the gap for larger intensities. 
This scenario is predicted correctly by the tight-binding model (1), however this model does 
not account for the existence of the second band, that defines the gap extent. On the other hand, 
the gap size limits the minimum width of self-trapped beams, see examples in Figs. 2(c,d). 
The width of the gap is smaller for weaker refractive index contrast and increases for larger 
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Fig. 2. (a,b) Bandgap spectrum of linear waves for index contrasts of 1.1 x 10 -4 and 
2.8 x 10 -4 , respectively. (c,d) Profiles of staggered gap solitons having the minimum width 
for cases (a,b), respectively. Shading marks index maxima, (e) Minimum width of the gap 
soliton vs. the refractive index contrast; (f) Efficiency of beam self-trapping calculated as 
the power fraction remaining in the 20 central waveguides at the output for an optimized 
input power (solid) compared with linear diffraction (dashed) vs. the refractive index con- 
trast. 



An m<3X . The calculated minimal width of the gap soliton, W = 3 J\x\\ \j/\ 2 dx/ J | \j/\ 2 dx, is plotted 
in Fig. 2(e) vs. the refractive index contrast. In the case of small refractive index contrast and 
a narrow band gap, the narrowest soliton spans over several waveguides. The situation changes 
when the contrast of index modulation increases, and the narrowest soliton is localized at a 
single waveguide. 

4. Crossover from self-defocusing to discrete self-trapping 

We study the crossover to discrete self-trapping by modeling the dynamics of an input Gaussian 
beam, which width is equal to the size of a single waveguide. In Fig. 2(f) we plot the relative 
power that remains in the central section of the array containing 20 waveguides (implemented 
with absorbing boundary conditions) after the propagation over a distance of many diffraction 
lengths (1000mm) vs. the index contrast of the lattice. We investigate the effect of nonlinearity 
(solid line) on beam dynamics with respect to linear diffraction (dashed line). We optimized the 
power of the input beam in order to maximize the power fraction which remains in the central 
section of the array. Our results demonstrate that for low index contrast the optimization of the 
input power for nonlinear propagation simply matches the linear limit, and the two curves co- 
incide in Fig. 2(f). In this regime, the nonlinear self-action results in increased beam spreading 
due to self-defocusing. As the refractive index increases, there appears a bifurcation from the 
linear regime. This corresponds to the formation of initially broad gap solitons, in agreement 
with the limitation on the soliton width shown in Fig. 2(e). As the minimum width of the soli- 
ton decreases and approaches that of a single guide for larger An max , self-trapping with almost 




Fig. 3. (a) Intensity distribution at the output facet of the array for linear propagation at low 
laser power (10 nW). (b) Corresponding intensity profiles: Solid - experimental measure- 
ment; shading - numerical solution of the full model [Eq. (2)]; crosses - lattice site am- 
plitudes calculated from the discrete model [Eq. (1)]. (c) Evolution of the beam intensity 
along the sample for a low input power simulated with Eq. (2). (d-f) Same as (a-c) for 
nonlinear propagation at high laser power (1 mW). 

100% efficiency becomes possible, as predicted by the discrete model (1). This feature repre- 
sents a transition from self-defocusing to discrete self-trapping as the index contrast exceeds a 
certain threshold. 

5. Experimental observation of single-site self-trapping with defocusing nonlinearity 

In previous experimental studies of soliton formation in permanent [5, 6] and optically- 
induced [8, 9] periodic structures with defocusing nonlinearity, the input excitation was spe- 
cially prepared to reflect the phase structure of the staggered solitons, where the neighboring 
sites are out-of-phase [see Figs. 2(c,d)]. This was achieved either by illuminating the structure 
with a single beam inclined at the Bragg angle [5, 8, 9] or by phase modulation of the initial 
beam [6] . However, observation of a soliton localization at a single waveguide was not reported. 

In order to study experimentally the discrete self-trapping at a single site supported by de- 
focusing nonlinearity, we designed an array of closely spaced optical waveguides produced by 
Titanium indiffusion into a mono-crystal lithium niobate wafer. The waveguide array has a pe- 
riod of 19/im and refractive index contrast An max = 2.8 xlO -4 , which was chosen above the 
threshold for a crossover to discrete self-trapping predicted in numerical simulations [Fig. 2(f)]. 
In the fabrication process, 100 A ofTi was deposited on the X-cut LiNb03 using electron beam 
evaporation. The Ti layer was then photolithographically patterned and etched in a buffered 
hydrofluoric acid solution. The diffusion was conducted at 1050°C for 3 hours in a wet oxygen 
environment. The waveguides were verified as single mode using a prism coupling technique. 
The array was then diced to a total length of 5cm and both facets were mechanically polished. 

The LiNbC>3 sample exhibits a strong photovoltaic effect which leads to the negative (self- 



Fig. 4. (a) Saturated camera image of the single channel localized state at the output of 
the array, zoomed at the dashed rectangle in Fig. 3(d). (b) Interferogram confirming the 
staggered phase structure of the output beam. 

defocusing) nonlinear response to laser beams at a photosensitive wavelength. In our experi- 
ments, we tightly focused an extraordinary polarized laser beam from a cw Nd: YVO4 laser into 
a single guide of the array by a microscope objective (x20). The input and output facets were 
monitored by two CCD cameras. The array was externally illuminated with white light in order 
to reduce the nonlinear response time of the photovoltaic material to less than a minute. At low 
laser power (~10 nW) the propagating beam experienced typical discrete diffraction [1], where 
at the array output most of the laser power was transferred into the neighboring waveguides and 
there was almost no light in the central guide [Fig. 3(a)]. The resulting intensity profile shown 
in Fig. 3(b, solid line) matches well the results of numerical simulations performed using the 
discrete model [Eq. (1)] (crosses) and the full model with periodic index modulation [Eq. (2)] 
(shading). The corresponding calculated propagation inside the array is depicted in Fig. 3(c). 

When the laser power is increased (1 mW) the defocusing photovoltaic nonlinearity leads 
to strong beam localization at a single waveguide [Fig. 3(d-Q], in a similar way as in focusing 
waveguide arrays [1], due to the universal nature of discrete self-trapping as discussed in Sec. 2. 

In order to confirm that the nonlinear state is indeed localized inside the Bragg reflection gap, 
and not in the total internal reflection gap, it is important to verify its staggered phase structure. 
For this purpose we allowed for saturation of the camera in order to detect the small but nonzero 
amount of light in the neighboring waveguides. The corresponding images are shown in Fig. 4 
where the two neighboring satellites are clearly visible. They are separated from the central 
waveguide by zero-intensity lines, which is an indication of their out-of-phase structure. To 
confirm the staggered phase of the localized state we interfere the output with an inclined broad 
reference beam. The corresponding interferogram is shown in Fig. 4(b), where a n shift of the 
interference fringes at the zero intensity lines is clearly observed. 

6. Conclusions 

We have studied the crossover between the beam self-defocusing and discrete self-trapping 
in waveguide arrays with defocusing nonlinearity. We have shown that an abrupt transition is 
observed in the beam dynamics with the increase of the refractive index. For a small index 
contrast, the beam experiences enhanced spatial spreading at higher input powers. However, 
when the contrast exceeds a critical value, beam self-trapping associated with the formation of 
staggered gap solitons becomes possible. We have demonstrated experimentally the generation, 
by single waveguide excitation, of strongly localized staggered states supported by defocusing 
nonlinearity in waveguide arrays, where the refractive index modulation was engineered to 
exceed the crossover threshold. 
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